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FOREWORD 


The  work  reported  here  was  oricinally  done  under  the  T-1S 
Subtask,  Comparison  of  Lifting  and  Ballistic  Reentry  Spacecraft 
for  Military  Purposes.  Results  of  the  initial  application  of 
the  derived  expressions  to  evaluate  glide  and  landing  character¬ 
istics  of  low  L/D  vehicles  are  described  in  IDA  Study  S-112 
(Aug  1963).  The  derived  linearized  expressions  have  been  in¬ 
corporated  in  the  RANGE  program  and  are  used  in  all  aircraft 
landing,  lifting-reentry,  and  horizontal- takeoff  launch  vehicle 
trajectory  calculations;  the  linear  approximations  are  funda¬ 
mental  to  all  other  RANGE  trajectories  involving  lift.  Subsequent 
IDA  documents  whose  results  make  use  of  these  routines  include 
P-237  (Jan  1966),  R-114  (Feb  1966),  P-427  (May  1968),  and  P-425 
(Jun  1969).  This  present  document  is  assembled  to  provide 
background  detail  for  these  aerodynamic  subroutines  which  are 
in  continuing  use  in  RANGE. 
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ABSTRACT 


The  dependences  of  lift  and  drag  coefficients  on  angle  of  attack 
for  bisymmetric  lifting  vehicles  are  synthesized  using  two  simpli¬ 
fying  assumptions: 

1.  the  axial  force  coefficient  is  independent  of  the  angle 
of  attack,  and 

2.  the  normal  force  coefficient  increases  linearly  with 
the  angle  of  attack. 

Good  agreement  with  experimental  data  is  found  for  the  value  of  the 
angle  of  attack  at  which  the  maximum  L/D  occurs  and  the  angle  of 
attack  for  maximum  CL  in  the  hypersonic  regime;  applicability  for 
subsonic  flight  is  limited  to  those  angles  of  attack  below  the  in¬ 
ception  of  flow  separation. 

Conditions  for  terminal  equilibrium  glide  and  landing  for  air¬ 
craft  are  derived  from  the  simplified  aerodynamic  characteristics. 

The  following  useful  rules  of  thumb  are  obtained  for  super  and 
hypersonic  vehicles  with  maximum  L/D  greater  than  1.0: 

1.  the  drag  at  maximum  L/D  is  very  closely  twice  the  zero-lift  drag, 

2.  the  L/D  at  maximum  CL  is  about  0.8, 

3.  the  angle  of  attack  for  maximum  CL  is  about  48  deg, 

4.  the  ratio  of  the  maximum  CL  to  the  CL  at  L/D  is 
approximated  by  0.5  +  L/Dmax>  and 

5.  the  ratio  of  the  normal  force  slope  to  the 
maximum  lift  is  closely  2.0  per  radian. 
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Reference  area 

Axial  force  coefficient 

Drag  coefficient 

Lift  coefficient 

Normal  force  coefficient 

Factor  relating  CN  and  its  a  dependence 

Vertical  force  coefficient  =  vertical  force/5gpVsA 

Drag  force 

Ratio  of  vertical  force  to  effective  weight 
Axial  force 

Flare  acceleration  in  g’s 

Lift  force 

Velocity 

Effective  weight  =  mass  times  local  acceleration  of 
gravity  less  the  centrifugal  force  due  to  velocity 

Angle  of  attack 

Initial  angle  of  attack 

Angle  of  attack  at  which  maximum  L/D  occurs 

Angle  of  attack  at  which  maximum  vertical  force  (or  CL 

Angle  of  attack  at  which  maximum  or  minimum  CN  occurs 

Flight  path  angle  with  horizontal,  positive  upwards 

Initial  or  desired  flight  path  angle 

Atmospheric  density 
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I.  INTRODUCTION 

Calculation  of  the  flight  path  of  a  point  mass  vehicle  acting 
under  the  influence  of  drag  and  lift  forces  requires,  in  addition  to 
the  more  commonly  available  dependences  of  these  forces  on  Mach  and 
Reynolds  numbers,  the  dependences  of  the  lift-drag  ratio  and  lift  or 
drag  on  the  angle  of  attack.  (With  a  point  mass  vehicle  the  drag  and 
lift  are  assumed  to  act  through  the  center  of  mass;  the  otherwise 
interesting  moments  are  assumed  to  have  been  cancelled  out  by  hypo¬ 
thetical  control  surfaces.)  A  description  is  given  here  of  the 
synthesis  of  these  aerodynamic  force  characteristics  for  a  vehicle 
configuration  simplified  but  sufficiently  general  to  generate  flight 
behavior  information  indicative  of  the  performance  of  most  real  vehicles 
in  the  angle  of  attack  range  0-20  deg  for  subsonic  flight  and  0-50  deg 
for  hypersonic  flight. 

The  subject  vehicle  is  required  to  have  two  basic  limitations. 

The  first  is  that  it  should  have  mirror  symmetry  in  the  two  perpendi¬ 
cular  planes,  vertical  and  horizontal,  whose  intersection  forms  the 
longitudinal  axis  of  the  vehicle.  Its  second  limitation,  really  an 
extension  of  the  symmetry  limitation,  is  that  stability  and  attitude 
control  are  to  be  supplied  by  control  surfaces  or  forces  giving 
negligible  drag  perturbation  or  mass  expenditure.  In  other  words  it 
should  be  a  bisymmetrical  body-wing  configuration  with  symmetry  not 


The  reproduction  of  this  Note  by  IDA  is  a  service  to  the  author 
and  the  professional  community  and  does  not  imply  IDA  endorsement 
of  the  views  expressed.  Its  release  outside  IDA  has  been  approved 
in  order  to  promote  discussion  of  the  ideas  presented. 
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only  left  to  right  but  also  with  symmetry  in  the  top  and  bottom  sur¬ 
faces  (airfoil)  and  having  a  method  of  attitude  control  that  does  not 
disturb  the  symmetry.  Bodies  of  revolution,  e.g.,  cylindrical  or 
conical  missile  bodies,  are  included  in  the  family  of  vehicles  covered 
by  these  limitations.  The  symmetry  limitation  is  to  assure  that  the 
force  normal  to  the  vehicle  is  zero  for  zero  angle  of  incidence  be¬ 
tween  the  longitudinal  axis  and  the  flight  path. 
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II .  BASIC  ASSUMPTIONS 


Two  assumptions  about  the  aerodynamic  force  characteristics  of  the 
subject  vehicle  are  required,  from  which  the  further  necessary  relations 
are  easily  derivable.  The  two  basic  assumptions  for  this  paper  are 
determined  empirically  from  inspection  of  a  wide  range  of  experimental 
results.  The  first  is  that  the  axial  force  coefficient,  the  relative 
aerodynamic  force  in  the  direction  of  the  vehicle’s  longitudinal  axis, 
is  a  constant  independent  of  the  angle  of  attack.  The  second  is  that 
the  normal  force  coefficient,  the  relative  aerodynamic  force  perpendicular 
to  the  longitudinal  axis,  has  a  quadratic  dependence  on  the  angle  of  at¬ 
tack.  In  equation  form  these  assumptions  are 

CA  =  constant  (1) 

CN  =  CN(X  o(l  -  a/2cO  (2) 

where  a,  is  the  angle  at  which  a  maximum  or  minimum  value  of  CN  occurs. 

The  first  assumption  is  good,  i.e.,  within  20%  in  the  first  20° 
of  a,  for  thick  delta  wings  (Ref.  1,  Mach  3-6),  excellent,  i.e., 
within  10%  for  a  <  20°,  for  cylindrical  bodies  (Ref.  2,  Mach  6.86), 
but  of  varying  applicability  for  sharp  cones  which  show  either  a  20% 
drop  in  CA  (Ref.  2,  Mach  2)  or  a  50%  increase  in  CA  (Ref.  5,  Mach  6.77) 
in  the  first  12°  of  a.  The  accuracy  of  the  assumption  becomes  poorer 
as  the  Mach  number  decreases  below  2,  but  remains  within  50%  at  a  <  12° 
even  at  low  subsonic  speeds  (Ref.  4). 


The  quadratic  CN  dependence  on  a  could  be  made  to  follow  the 
data  within  a  few  percent  up  to  a  »  45°  throughout  the  whole  range  of 
Mach  numbers  if  a,  were  allowed  to  vary  with  M.  Generally  a,  would  have 
a  large  negative  value  leading  to  a  CN  vs  a  curve  that  is  slightly  con¬ 
cave  upward  (Refs.  1,  2,  and  5).  The  value  of  a,  would  become  closer 
to  zero  as  the  Mach  number  increases,  leading  to  a  more  nonlinear  CN 
vs  a  curve. 

A  linear  approximation  to  assumption  2  (a,  =  ±  ®)  is  generally 
acceptable  (deviation  less  than  10%)  up  to  a  =  20°  and  sometimes  higher. 
Newtonian  impact  theory  predicts  an  increase  in  CD  with  a  (Ref.  10), 

ac„ 

*r- 5C<- 

that  is  as  much  as  30  to  40%  greater  than  supported  by  experiment 
(Refs.  11  and  12)  and  about  50%  greater  than  that  given  by  applying 
the  assumption  CA  =  constant  to  the  differential  of  equation  4  below. 
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III.  DERIVATION  OF  RELATIONS 

The  basic  relationship  is  the  transformation  from  the  body  re¬ 
ference  system  of  forces,  i.e.,  the  axial  and  normal  forces  on  the 
vehicle,  to  the  flight-path  reference  system  defined  by  the  lift  force 
perpendicular  to  the  flight  path  and  the  drag  force  parallel  with  the 
flight  path.  From  the  force  diagram  in  Figure  1 


CL 

~ 

CN  cos  a  - 

CA  sin  a 

(3) 

Co 

= 

CN  sin  a  + 

Ca  cos  a 

(4) 

or  CL/C. 

= 

tan  (tan-1 

Cn/Ca  -  a) 

(5) 

transformation  is 

Cn 

= 

CL  cos  a  + 

C0  sin  a 

(6) 

Ca 

= 

C0  cos  a  - 

CL  sin  a 

(7) 

or  CN/CA 

= 

tan  (tan"1 

CL/Co  +  a) 

(8) 

Differentiation  of  equation  5  with  respect  to  a  (with  substitution  of 
equation  2)  and  setting  the  result  equal  to  zero  gives  a  preliminary 
relation  for  ax  at  which  maximum  L/D  occurs 


ai  = 


/( CNa/CA)  (1  -  ax/a.)  -  1 
(CNa/CA)  (1  -  ax /2a.) 
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Vertical 


FIGURE  1 .  Force  Diagram  for  a  Gliding  Aircraft  Showing  Relationship  of  Flight-Path- 
Oriented  Force  Coefficients,  C^and  C^,  with  Body-Oriented  Force 

Coefficients,  C^  and  C^,  and  with  an  Earth-Oriented  Force  Coefficient,  Cy 
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or 


(CNa/CA)  (1  -  c^/a.)  -  1  =  (CN/CA)2 


a  -  ax 


(10) 


A  most  useful  defining  parameter  of  the  vehicle  is  the  maximum  L/D 
value.  Substituting  equation  8  in  equation  10,  we  obtain  as  an  inter¬ 
mediate  result  in  terms  of  L/D 

max 

sec2  (tan"1  L/Dmav  +  ax)  (11) 


Cng/C*  "  (i  .  cli/O 


which,  when  substituted  in  equation  9,  gives 

tan  (tan-1  L/D _  +  ax)  (1  -  ai/aB) 


ai  - 


max 


(12) 


sec2  (tan-1  L/Dmax  +  o.i)  (1  -  ai/2a.) 


which  can  be  finally  simplified  to 

ai  =  %  sin  2  (tan"1  L/Dmav  +  ax)  (1  -  ai/a,) 

(1  -  ai/2a.) 


max 


(13) 


Using  this  value  of  ax  and  a  combination  of  equations  2  and  8,  a 
slightly  simpler  relation  than  equation  11  is  obtained 


p  /r  _  ton  (tan"1  L/Dmax  +  cQ 
Na/C*  af  (1  -  ai/2a,) 


(14) 


Differentiating  equation  3  with  respect  to  a  (with  substitution  of 
eqs.  1  &  2)  and  setting  equal  to  zero  gives  a  relation  for  aa  at  which 
maximum  CL  occurs 

Cn„/Ca  (1  -  Og/o.)  -  1 

(15) 


da  tan  da  = 


CNa/CA  (  1  -  Oa/20.) 
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For  glide  at  an  angle  y0  negative  below  the  horizontal,  the  vertical 
force  coefficient  (in  an  earth-oriented  system  rather  than  body-oriented 
or  f light -path-oriented ,  Fig.  1)  is 

Cv  --  CN  cos  (Yo  +  a.)  -  CA  sin  (y0  +  a)  (16) 


Maximizing  this  with  respect  to  a  gives  an  alternative  equation  for 
aB,  of  which  15  is  a  special  case  for  y0  =  0, 


a,  tan  (y0  +  Oa)  = 


CnqA-a  (1  -  Qfl/a.,)  -  1 

CNa/CA  (1  -  09/20.) 


(17) 


This  gives  the  angle  of  attack  at  which  the  maximum  vertical  force  is 
available  to  act  against  the  vehicle  weight. 


One  additional  relationship  is  obtained  from  inspection 
of  equation  8.  If  tan”1  (L/D  +  a)  is  greater  than  tt/2,  then  CN/CA  is 
negative,  which  is  illogical.  Consequently,  we  have  a  constraint  on 
the  logical  values  of  a  for  a  given  value  of  L/D, 

a  £  tt/2  -  tan”1  L/D  (18) 


The  region  of  a  greater  than  this  quantity  is  an  excluded  region. 

Figure  2  shows  a  set  of  plots  of  relations  5,  13,  and  18  describing 
the  dependences  of  L/D  on  a  for  vehicles  with  a,  =  »  and  different 
values  of  L/Dm3v ,  the  locus  of  the  maxima  of  L/D  vs  a  curves  for 
vehicles  with  a,  =  °°,  and  finally  the  boundary  of  the  excluded  region. 
Figure  3  gives  the  equivalent  dependences  to  L/D  on  Q.  from  equations 
3  and  15,  with  CL  defined  in  terms  of  the  ratio  to  Q.max. 
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IV.  EQUILIBRIUM  GLIDE 


«r 


A  particular  vehicle  requires  the  value  of  one  more  parameter  in 
addition  to  L/D,^  to  define  completely  its  aerodynamic  force  character¬ 
istics  for  a  point-mass  trajectory  calculation.  That  parameter  will 
most  likely  be  its  axial  force  coefficient,  i.e.,  its  zero-lift  drag 
coefficient,  for  subsonic  applications  or  its  maximum  lift  coefficient 
for  hypersonic  applications.  For  a  vehicle  in  subsonic  glide,  the 
value  of  Ck  can  be  obtained,  for  example,  from  the  vehicle’s  terminal 
velocity  at  a  constant  glide  angle  (equilibrium  glide  condition),  in 
the  following  manner: 

For  an  equilibrium  glide  at  an  angle  y0  which  is  measured  from  the 
horizontal  and  negative  downward,  the  resultant  of  the  aerodynamic  forces 
in  Fig.  1  is  equal  and  directly  opposite  to  the  effective  weight  W tff,  i.e., 


L 

=  VJ,tt  cos  Yo 

(19) 

D 

=  -W,  f  t  sin  Yo 

(20) 

or 

L/D 

=  -cot  Yo 

(21) 

This  relationship  is  the  same  as  the  equality  in  18  with  a  replaced  by 

-  Y0 •  With  L/D  in  place  of  L/D,  this  expression  also  gives  the 
rnsx 

minimum  equilibrium  glide  angle.  The  boundary  of  the  excluded  region 
in  Figure  1  therefore  gives  the  L/D  required  for  equilibrium  glide  at 
an  angle  y0  reading  -  Yo  as  the  abscissa  instead  of  a,  or  the  minimum 
glide  angle,  Ymin>  for  vehicles  with  different  L/Dmax- 
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From  equation  21  and  equation  8  we  can  derive  the  angle  of  attack 


of  the  vehicle  in  this  equilibrium  glide 
ao  =  -Yo  -cot"1  kCN  /CA)a0  (  1 


-It.  >] 


(22) 


where  the  solution  a0  <  a*  is  chosen.  Finally  from  equation  7  and 
equations  19  and  20 


C*  = 


Co  cos 


a0  -  CL  sin  ae 

=  -  wjX1.  fin  Y>  cos  o.  -  W.M  Sin  a. 

1/2  pV2A  1/2  pV3A 

=l72^V5ASln(-Y° 


(23) 


If  L/Dmax  and  CA  are  given,  then  the  equilibrium  glide  velocity 
is  found  from  equation  23, 


1/2 


V.  RECAPITULATION  FOR  LINEAR  DEPENDENCE  OF  CN  ON  a 

In  subsonic  glide  situations,  only  angles  of  attack  less  than 
20°  are  of  general  interest.  In  this  region  of  values  of  a  a  linear 
dependence  of  CN  on  a  may  be  used  with  little  error,  giving  some  simpli¬ 
fication  to  the  equations.  The  simplified  equations  are  listed  below 
in  the  order  of  their  calculation  with  a  description  of  their  method  of 
solution.  The  corresponding  FORTRAN  subroutines  are  listed  in  Appendix  A. 

The  angle  of  attack  for  L/D  „ 

3  max 

\  sin  2  (tan"1  L/Dmax  +  ax)  (13*) 

is  a  transcendental  equation.  It  is  solved  by  iteration  from  a  first 
guess  that  is  very  close 

■  I  <"'2  -  tan“  ■  -  I  Vin  <24> 

Each  iteration  input  after  the  first  is  the  average  of  the  input  and 
output  of  the  previous  iteration. 

The  ratio  of  the  normal-force-coefficient  slope  to  the  axial -force 
coefficient  is  a  simple  direct  equation 

cN  /cA  =  tan  L/Dmax  1  (14») 

a  ax 
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The  angle  of  attack  for  maximum  vertical  force  (maximum  CL  for 


horizontal  flight) 


cNa/cA  -  1 

(CNc/C*')  tan  (y0  +  a*) 


(17’) 


is  another  transcendental  equation  to  be  solved  by  iteration.  The  first 
guess  is 


a*  =  0.86  -  0.64  y0  (25) 


Each  iteration  input  after  the  first  is  a  value  only  one-quarter  of  the 
way  from  the  input  to  the  output  of  the  previous  iteration,  to  guarantee 
convergence  for  an  otherwise  divergent  process. 

With  the  above  equations  as  programmed  in  the  subroutines  in 

Appendix  A,  values  are  calculated  for  ax ,  Og,  CN  /CA  ,  CL  /CA  ,  CL  /CA  ; 

ol  ^  max 

CLmax/cD#>  and  Ymin>  and  are  listed  in  Table  1  for  the  interesting 

range  of  values  of  L/Dmax.  The  ratios  CLmax/CLi ,  CN(i/CLmax,  CWCA , 

and  C.  /CD  are  plotted  in  Fig.  4.  The  values  in  the  Figure  are 
max  2 

more  limited  in  their  range  than  the  values  in  the  table,  and  indeed 

the  latter  three  of  them  are  essentially  constant  for  L/D  ^  values 

max 

greater  than  1.0. 

With  the  linear  CN  dependence,  the  angle  of  attack  for  equilibrium 
glide  becomes 

a0  =  -  Yo  -  tt/2  +  tan"1  ((CNa/CA)  ae)  (22’) 

The  solution  to  this  equation  c_n  best  be  shown  graphically  with  the 
aid  of  Figure  2.  Reading  up  from  -  Yo>  as  the  abscissa,  to  the  boundary 
curve  gives  the  required  L/D.  Reading  across  with  this  L/D  there  are 
in  general  two  intersections,  or  none,  with  an  L/D  vs  a  dashed  curve. 
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TABLE  I 


Aerodynamic  Force  Characteristics  of  Generalized 
Vehicles  with  Different  L/Dmax 


l/dmv 

ax 

a* 

.100 

25.307 

28.274 

1.361 

.200 

28. 101 

34.108 

1.675 

.300 

28.643 

37.  621 

2.024 

.400 

28.209 

40.03S 

2.421 

.  500 

27.289 

41.796 

2.874 

.600 

26.132 

43. 124 

3.389 

.700 

24.875 

44.150 

3.968 

.800 

23.602 

44.955 

4.616 

.900 

22.359 

45.  597 

5.336 

1.000 

21. 1 73 

46.116 

6.128 

1.100 

20.056 

46. 538 

6.994 

1  .200 

19.014 

46.88  5 

7.936 

1.300 

18.046 

47.175 

8.955 

1.400 

1  7. 1  50 

47.417 

10.050 

1.500 

16.323 

47  .  622 

11.223 

1.600 

15.558 

47.796  • 

12.475 

1.700 

1  4. 8  52 

47.945 

13.804 

1.800 

14. 199 

48.074 

15.213 

1.900 

13.594 

48.186 

1  6.  700 

2.000 

13.034 

48.283 

18.267 

2.100 

12.  513 

48.369 

19.912 

2.200 

12.029 

48.444 

21.637 

2.300 

1 1 . 579 

48. 51 1 

2$. 441 

2.400 

1 1 . 1  58 

48. 570 

.  25.325 

2.500 

10.765 

48  .  623 

27.289 

2.  600 

10.397 

48.671 

29.332 

2.  700 

10.052 

48.714 

31.455 

2.S00 

9.728 

48.752 

33.657 

2.900 

9.423 

48.787 

35.940 

3.000 

9.136 

48.819 

38.302 

3. 100 

8.865  ' 

48.847 

40.744 

3.200 

8.609  . 

48.874 

43.265 

3.300 

8.367 

48.898 

45.867 

3.400 

8.138 

48.920 

48. 548 

3.  500 

7.920 

48.940 

51.310 

3.600 

7.714 

48.959 

54.151 

3.  700 

7.517 

48.976 

57.072 

3.800 

7.331 

48.992 

60.073 

3.900 

7.152 

49.007 

63.154 

4.000  ' 

6.983 

49.021 

66.315 

4.  5(10 

6.239 

49.077 

83.319 

5.000 

'  5.636 

49.117 

102.322 

5.500 

5.138 

49.147 

123.324 

6.000- 

4.720 

49.170 

146.325 

6.500 

4.365 

49.188 

171.326 

7.000 

4.0  58 

,  49*203 

198.327 

7.500 

3.792 

49*21 5 

227.328 

8.000 

3.558 

49.224 

258.329 

b.  500 

3.351 

49.232 

291.329 

9.000 

3.167 

49.239 

326.330 

9.5(10 

3.002 

49.244 

363*  330 

10.000 

2.853 

49.249 

403*331 

cLi/cA 

1  ~ >mlr 

.116 

■  .ii e ‘ 

•  09o 

.254 

.  .265 

.  1  9  i 

.409 

.442 

.276 

.  578 

.  652 

.352 

.  758 

.897 

.418 

63.43 

.947 

1.178 

.476 

1.142 

1.497 

.  52  6 

1.342 

1.857 

.  568 

1 . 545 

2.257 

.  6»<  4 

1.750 

2.698 

.  63  5 

45.00 

1.957 

3.182 

.  661 

2.164 

3.709 

.  684 

2.372 

4.278 

.703 

2.  580 

4.892 

.719 

2.  787 

5.549 

.  733 

33.69 

2.995 

6.2  50 

.746 

* 

3.202 

6.995 

.756 

3.409 

7.78  5 

.766 

3.  61  6 

8.619 

.  774 

3.823 

9.497 

.781 

26.55 

4.029 

10.420 

.  788 

4.235 

11.387 

.793 

4.440 

12.399 

.  798 

4.  645 

13.456 

.80  3 

4.8  50 

14.557 

.807 

21.82 

5*055 

15.704 

.81  1 

5.259 

16.894 

.814 

5.464 

18.130 

.817 

S.  667 

19.410 

.  820 

5.871 

20.736 

.822 

18.43 

6.075 

22.106 

.824 

6.278 

23.520 

.827 

6.481 

24.980 

.628 

6.  685 

26.484 

.830 

6.887 

28.034 

.832 

15.94 

7.090 

29.628 

.633 

7.293 

31.267 

.835 

7.495 

32.951 

.836 

7.698 

34. 679 

.637 

7.900 

36.453 

.838 

14.04 

8.91 1 

45.993 

.843 

12.52 

9.919 

56.655 

.846 

11.31 

10.926 

68.439 

.845 

10.30 

1 1.932 

81.345 

.  h  SI’ 

9.46 

12.937 

95.373 

.  8  52 

8.75 

13.941 

110.523 

.  8  So 

8.13 

14.945 

126.795 

.854 

7.59 

15.948 

144. 190 

.  8  5  j 

7.13 

1  6.951 

162.706 

.855 

6.71 

1  7.954 

182.345 

.856 

6.34 

18.956 

203.105 

.856 

6.01 

19.959 

224. 9B8 

.8  57 

5.70 

Ratio  of  coefficients 


« 


The  left,  lower  a,  intersection  is  of  greater  interest  because  it  re¬ 
presents  an  aircraft  with  a  greater  velocity  excess  above  stall.  Un¬ 
fortunately,  a  simple  iterative  process  converges  on  the  other  inter¬ 
section.  Consequently,  equation  22'  is  solved  by  trial  and  linear 
interpolation  with  the  first  two  trial  values  1/2  ax  and  1/4  ax . 

The  axial  force  coefficient,  or  zero-lift  drag  coefficient,  is 
given  by  equation  23  without  change.  Since  the  atmospheric  density 
enters  in  equation  23,  the  above  set  of  equations  for  equilibrium  glide 
will  give  different  values  at  different  altitudes. 

In  subsequent  flight  to  an  altitude  with  a  different  atmospheric 
density,  the  vehicle  will  require  an  angle  of  attack  to  maintain  zero 
vertical  acceleration  (which  is  not  necessarily  the  same  as  equilibrium 
glide  as  defined  here  by  equations  19  and  20)  given  by  solution  of  the 
following  equation  by  trial  and  linear  interpolation 

(CNa/CA)  a  F*  cos  (y  +  a)  -  FA  sin  (y  +  a)  -  W,f f  =  0  (26) 

The  fir3t  guess  for  a  is 

a  =  W.ff/((CMa/CA)  Fa)  (27) 

If  a  given  glide  slope  Yo  is  to  be  maintained  instead  of  zero 
vertical  acceleration,  the  effective  weight,  W, ff,  is  replaced  in 
equations  26  and  27  by  fW, ff  where  f  has  a  value  near  1  given  by 

f  ~  1  -  (y-Yo)  (28) 

For  a  flare  maneuver  requiring  a  pullup  acceleration  9r^are>  the 
expression 

f  =  1  +  »fUre  (29) 

gives  very  closely  the  proper  flare  acceleration. 
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For  horizontal  flight  after  a  flare,  equation  28  is  used  with 

Ye  =  o. 

Using  the  above  relations,  the  performance  in  the  glide,  flare,  and 
landing  regimes  has  been  calculated  with  the  IDA  program  RANGE  for  general 
bisymmetric  lifting  vehicles  with  different  wing  loading  and  L/Dmax.  The 
results  have  been  reported  in  Ref.  6. 
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VI.  COMPARISON  OF  LINEARIZED  AERODYNAMIC 
CHARACTERISTICS  WITH  EXPERIMENTAL  DATA 

Wind-tunnel  tests  of  some  bisymmetric  thick  delta-wing  models  at 

Mach  3,  4.5,  and  6,  are  reported  in  Ref  1.  Equation  13’  predicts  the 

angle  of  attack  a*  giving  L/D  to  about  1°  accuracy  out  of  15°. 

max 

Equations  17’  and  14’  predict  the  CLmax  angle  of  attack  a®  about  10 
percent  low  for  the  most  non-linear  CN  vs  a  curves. 

The  X-20  Dyna-Soar  vehicle  could  hardly  be  called  symmetrical 
above  and  below  the  horizontal  plane  but  it  is  still  interesting  to 
compare  the  calculations  (linearized  equations  13',  14',  and  17')  with 
its  behavior  measured  in  wind-tunnel  tests  (Refs.  4  and  7).  The  table 
below  gives  the  comparison.  The  angles  of  attack  are  measured  above 
the  zero-lift  attitude. 

Subsonic  (L/Dm_v  =  4)  Hypersonic  (L/Dn=v  =  2) 

max  illaX 


2± 

2± 

2 L 

da 

Measured 

8° 

O 

O 

12° 

46' 

Calculated 

/ 

49° 

13° 

48' 

The  relatively  poor  subsonic  agreement  in  Og  is  a  result  of  flow  sepa¬ 
ration. 

Lift/drag  data  are  plotted  in  Fig.  5  for  the  above  vehicles  and  for 
three  others  superimposed  on  the  linearized  L/D  vs  a  curves  from  Fig.  2. 
The  curves  reproduce  the  dependences  shown  by  these  data  generally  within 
2  deg  of  angle  of  attack. 
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FIGURE  5.  Comparison  of  Linearized  Aerodynamic  Characteristics  with  Experimental  Data 


VII.  CONCLUSIONS 


« 


Useful  relations  for  estimating  the  aerodynamic  force  character¬ 
istics  of  general  bisymmetric  lifting  vehicles  in  all  regimes  of 
flight  can  be  derived  from  the  assumptions  of  a  linear  variation  of 
the  normal  force  coefficient  with  angle  of  attack  and  an  independence 
of  the  axial  force  coefficient  on  the  angle  of  attack.  Good  agreement 
with  wind-tunnel  measurements  is  observed. 

Some  extremely  simple  rules  of  thumb  are  found  from  the  linearized 
aerodynamic  characteristics  relations.  In  the  range  of  maximum  L/D 
above  1.0  the  following  approximate  rules  hold: 

1.  The  drag  at  the  angle  of  attack  for  maximum  L/D  is  very 
nearly  twice  the  zero-lift  drag,  or 

CD  /CA  «  2 

l 

and 

CS/C.  „  2  L/Dmax 

2.  The  lift-drag  ratio  at  maximum  lift  coefficient  is  about 
0.8  (cf.  Ref.  8),  or  L/D  (?  CLmax  «  0.8 

3.  The  angle  of  attack  for  maximum  lift  coefficient  is  about 
48  deg,  or 

as  48  deg 
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4.  The  ratio  of  the  maximum  lift  coefficient  to  the  lift  co¬ 


efficient  at  maximum  L/D  is  given  by 


^max/S  ~0*5  *  L/Dmax 


and  from  the  first  rule  above 


W2*  “  L^Dmax  +  2  (L/Dmax> 

5.  The  ratio  of  the  normal  force  slope  to  the  maximum  lift  de¬ 
creases  only  from  2.25  per  radian  at  L/D  =  1.0  to  1.80  per 

max 

radian  at  L/D  =  5.0.  so 
max 


C»c/CW  -  2  rad 


-1 


For  lifting  bodies,  the  hypersonic  maximum  lift  coefficient  is 
about  0.6  (Refs.  1,  7,  and  8),  referenced  to  the  planform  area, 
so 

CN  «  1.2  rad"1  for  lifting  bodies. 

For  cones,  the  hypersonic  maximum  lift  coefficient  is  about 
1.0  (Refs.  2  and  5),  referenced  to  the  base  area,  so 

CN„  «  2.0  rad”1  for  cones. 

UL 
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APPENDIX  A 


Calculation  of  Aerodynamic  Force  Characteristics  Using  Equations  13’, 
14’,  24',  17*,  and  25*;  and  Program  to  Generate  Table  I. 

ivAOAN  G  =  .0174532925 
110  DO  7  LDMAX  «  1  #40 
.  1  ALPH1  =  0 

. >  CMACA  =  0 

1  FLDMAX  =  LDMAX 

i.-ii-i  FLDMAX  =  0.1  *  FLDMAX 

call  flt<fldmax#  alphi#  cnaca) 

1;.,..  ALPH2  =  ALP (CNACA# ZERO) 

16fc  CL1CA  =  CNACA  *  ALPH1  *  CO SF( ALPHI)  -  SINFCALPH1) 

i  Via  CLMCA  ■  CNACA  *  ALPH2  *  C0SF(ALPH2>  -  SINFCALPH2) 

1  V/.  CLMCD  =  TANF CATANF C CNACA  *  ALPH2)  *  ALPH2) 

.  V .i  ALPH1D  =  ALPH1/RADAN6 

i iwfi  AL1-H2D  =  ALPH2/RADANG 

210  7  PRINT  9#  FLDMAX,  ALPHI D#  ALPH2D#  CNACA#  CL1CA#  CLMCA#  CLMCD 

226  9  FORMAT  C7F10.3) 

230  END 


(13*) 


(14’) 


(24’) 


function  alp<x#y> 

4L3  PI  =  3.1415926536 
-v.  IF  (  A  )  6#  6#  7 

-.35  1  SY  ■  E  ♦  Y 

4*,.v  *1  ■  C  X  -  1  .  )/(X  *TANPC ZY ) )-- . 

45.5  IF  (ABSF(Z1  -  Z)  *  1  •  E  -  4)  3#  3#  2 
460  2  2  =  Z  ♦  (Z1  -  Z>  *  0*25 
470  GOTO  1 

4iv./  3  2  *  Z  ♦  (21  -Z>  *  0.25 
490  IF  (Z  -  PI/2.)  5#  5#  4 
5 ;>k>  4  c  «  PI/2. 

5  ALP  =  Z 

523  IF  CABSFCY)  -  .01)  9#  9#  8 

530  6  Z  »  0.86  -  Y  *  0.64 - r~ . — 

64w  GO  TO  1 

545  7  2  =  0.86  -  Y  *A 

.>*,6  GO  TO  1 

u50  8  A  =  CZ  -  0.86)/(-Y> 

560  RETURN 
„70  9  A  »  0.64 
575  END 


(17’) 


(25’) 


240  SUBROUTINE  FLT<X,Y#C) 
w 5U  PI  a  3.1415926536 

is  wi#  c  =  ATANF  (X) 

70  IF  CB)  4#  4#  5 

1  Y 1  =  0.5  *  SINFC2.  *  (Z  ♦  Y>)  — 
ivu  IF  (AbSF(Yl  -  Y)  -  l.E-4)  3#3#2 

2  Y  =  < Y 1  ♦  Y)/2. 

.  fj  go  to  i 

32.,  3  V  =  < Y 1  ♦  Y )/2. 

c  =  Tanfcz  ♦  Y)/Y . 

24..  b  =  Y/CPI/2.  -  Z) 

RETURN 

:.;G  4  Y  =  (PI/2.  -  Z)  *  0.48897 - 

60  TO  1 

.  5  Y  =  (PI/2.  -  H)  *  B 

C  *  .7  GO  TO  1 
END 
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The  dependences  of  lift  and  drag  coefficients  on  angle  of  at¬ 
tack  for  bisymmetric  lifting  vehicles  are  synthesized  using  two 
simplifying  assumptions:  (1)  the  axial  force  coefficient  is  inde¬ 
pendent  of  the  angle  of  attack,  and  (2)  the  normal  force  coefficient 
increases  linearly  with  the  angle  of  attack.  Good  agreement  with 
experimental  data  is  found  for  the  value  of  the  angle  of  attack  at 
which  the  maximum  L/D  occurs  and  the  angle  of  attack  for  maximum  Cu 
in  the  hypersonic  regime;  applicability  for  subsonic  flight  is  lim¬ 
ited  to  those  angles  of  attack  below  the  inception  of  flow  separa¬ 
tion. 

Conditions  for  terminal  equilibrium  glide  and  landing  for  air¬ 
craft  are  derived  from  the  simplified  aerodynamic  characteristics. 

With  the  linearized  normal  force,  the  following  useful  rules  of 
thumb  are  obtained  for  super  and  hypersonic  vehicles  with  maximum 
L/D  greater  than  1.0:  (1)  the  drag  at  maximum  L/D  is  very  closely 

twice  the  zero-lift  drag,  (2)  the  L/D  at  maximum  CLis  about  0.8,  (3) 
the  angle  of  attack  for  maximum  CLis  about  48  deg,  (4)  the  ratio  of 
the  maximum  CLto  the  CLat  L/D  max  is  approximated  by  0.5  +  L/Dmax> 

and  (5)  the  ratio  of  the  normal  force  slope  to  the  maximum  lift  is 
closely  2.0  per  radian. _ 
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